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ABSTRACT 

We present 15 - 20 /im long-slit spectra, from the Infrared Spectrograph (IRS) 
on Spitzer, of NGC 7023. We observe recently-discovered interstellar emission 
features, at 15.9, 16.4, 17.0, 17.4, 17.8, and 18.9 /im, throughout the reflection 
nebula. The 16.4 /im emission feature peaks near the photodissociation front 
northwest of the star, as do the aromatic emission features (AEFs) at 3.3, 6.2 
and 11.3 /im. The 16.4 /im emission feature is thus likely related to the AEFs 
and radiates by non-equilibrium emission. The new 18.9 /im emission feature, 
by contrast, decreases monotonically with stellar distance. We consider candi- 
date species for the 18.9 /im feature, including polycyclic aromatic hydrocarbons, 
fullerenes, and diamonds. We describe future laboratory and observational re- 
search needed to identify the 18.9 /im feature carrier. 

Subject headings: dust, extinction — infrared: ISM — ISM: individual (NGC 
7023) — ISM: lines and bands — ISM: molecules — reflection nebulae 



1. Introduction 

The spectrum of aromatic emission features (AEFs) between 3.3 and 12.7 /im has been 
previously well characterized by ground-based, airborne, and space-based spectroscopy. The 
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5-2 /im region of AEF sources, however, has been far less examined. iBeintema et al. 



(119961 ) obtained 2.4 - 45 /im Infrared Space Observatory (ISO) spectra of three bright 



AEF sources, using the Short Wavelength Spectrograph (SWS). Th ey found new emissio n 



features in the 15 - 20 /im regio n at 17.0 (b r oad) and 17.9 /mi. iMoutou et al.l (120001 ). 
Van Kerckhoven et al.l (120001 ) , and ISturm et al.l (120001 ) discovered a new emission feature at 
16.4 jum in a variety of sources, with ISO-SWS spectra. The 16.4 /mi feature was iden- 



(Moutou et al. 


2000) or the C- 


gs ( 


Van Kerckhoven et al. 


2000) 



covered a new 15.8 /m i emis sion feature in ISO-SWS spectra of five star-forming galaxies. 



Van Kerckhoven et al.l (120001 ) detected a broad plateau of emission at 15 - 20 jum in sever al 



sources, corresponding to the broad 17.0 /im feature discovered by IBeintema et ajj (119961 ). 



The Infrared Spectrograph (IRS; lHouck et al.ll2004l ) on Spitzer (jWerner et al.ll2004al ) has 
the sensitivity to detect si gnificant spect r al subs tructure in th e 15 - 20 /im region not evident 



in previous observations. I Werner et al.l (j2004bl ) (hereafter iPaper ll ) have discovered new 



spectral features at 17.4 and 19.0 /mi, in addition to previously detected narrow features at 
15.8, 16.4 and 17.8 /mi and a broad feature at 17.0 /mi, in Spitzer IRS spectra of NGC 7023. 
We describe in §3.3 how we use higher resolution spectroscopy to improve the wavelengths 
for these features; as a result, the '15.8' and '19.0' /im features become the '15.9' and '18.9' 
/tm features. These spectral features at 15.9, 16.4, broad 17.0, 17.4, 17.8, and 18.9 /mi have 



(Armus et al. 


2004 




Brandl et al. 


2004 




Morris ct al. 


200^ 


1; 


Smith et al. 


2004 




Hiedon et al. 


20061; 


Eneelbracht et al. 


20061; 


TaDDe!l2006: 


Panuzzo et al. 


2006: 


Smith et al. 


120061). including 


the new 18.9 /im feature which has been detected in several g 


alaxies JPanuzzo et al. 


2006; 



Smith et al.l 120061 ) . 



We extend here the work of IPaper ll . by presenting additional Spitzer IRS spec- 
troscopy of NGC 7023, a well-studied photodissociation region (PDR) and reflection neb- 
ula. NGC 7023 is illuminat ed by the Herbig Be star HD 200775, at a distance of 430 pc 
(Ivan den Ancker et al .11 1 9 9 71 ) . This extension of our previous work allows us to better char- 
acterize the spatial distributions and emission mechanisms of the features in the 15 - 20 /im 
region. 



2. Observations 



We describe our observational technique in lPaper ll . We improve the data reduction and 
examine here, in further detail, the Short-High (SH) IRS spectrum (9.9 - 19.6 /im; R = A/AA 
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= 600) of NGC 7023 at one pointing position previously presented in iPaper ll . obtained on 
2003 September 19. The SH integration time was 6 s, repeated twice. We extract the SH 
spectrum from the entire aperture (4'/7 x 11'.'3). We reduce three new slit pointing positions 
(24 additional spectra in NGC 7023) of second-order Long-Low (LL2) long-slit IRS spectra 
(14.0 - 21.3 //m; R = 80 - 128), obtained on 2003 October 1. The LL2 integration time was 
30s, repeated 3 times. We re fer to the se new slit pointing positions as slits A, C, and D. In 
this notation, the LL2 slit in IPaper ll is slit pointing position B. 



Figure [T] illustrates the location of slit pointing positions A, B, C, and D. We extract 
eight spectra from each slit, at evenly-spaced spatial positions we refer to as positions 1 
through 8, labeled from southeast to northwest (see Fig. [I]). All of the slits A, B, C, and 
D were observed at a slit position angle 306.4° east of north. Slit B crosses HD 200775 at 
position 3 (spectrum L3 in lPaper ll ). Slits A and C are shifted relative to B by 6'.'4 southwest 
and 19'.' 2 northeast, respectively, in a direction perpendicular to the slit. Slit D is shifted 
southeast by 45", relative to C, in a direction parallel to the slit. We extract spectra from 
fixed pixel ranges (extraction box of 2 x 3 pixels, or 10 / /2 x 15'/3) along the 2 pixel wide slit 
of LL2, producing eight spectra for each slit pointing position. We re-reduce Slit B, both 
for consistency with the other spectra presented in this paper and because of improvements 
in the calibration; this changes the intensities of Slit B spectra by ~7%. We present here 
an analysis of a total of 32 spectra. All of our data was obtained prior to nominal Spitzer 
operations (pre-nominal status), so our intensity, wavelength, and pointing uncertainties are 
larger than for nominal (current) Spitzer operations. The SH and LL2 absolute positions are 
thus uncertain by ±5". 

We reduce the data a nd extract spectra u sing Cornell IRS Spectroscopy Modeling Anal- 
ysis and Reduction Tool (IHigdon et al.l 120041 ) . All data are hand-calibrated by applying the 
spectrum and model template of point source calibrator HR 8585 (a Lac) to the intermediate 
un-flat-fielded product of the SSC pipeline. We subtract blank sky spectra only from the 
SH spectrum, observed as Spitzer cooled, to remove thermal emission from the then warm 
(~45 K) baffles of the telescope. We apply the current extended source correction to the 
LL2 spectra. The LL2 spectra have an absolute uncertainty in their intensities of ~15%, 
due to their being obtained very early in the Spitzer mission, before the detector bias was 
changed. 

Figure [1] shows that three of the four LL2 slits we observe intersect the filaments 
seen to the northwes t of HD 200775. The filaments, seen in Figure HJ trace a PDR front 
(ILemaire et al.lll996l ). All four LL2 slit s also co ver low-intensity regions to the southeast of 
HD 200775. The LL2 slit illustrated in lPaper ll is shown in Figure [Ha s Slit B, and it inter- 
sects both the northwest filaments and HD 200775. The SH spectrum (IPaper ll ) is obtained 
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0"W 29"N of HD 200775, between the star and the PDR front. 



3. Results 
3.1. Observed Features at High Resolution 

Figure El illustrates our SH spectr um (15.1 - 19.5 /un; R = 600) in NGC 7023. This 



is a part of the SH spectrum shown in iPaper ll . although we improve the reduction, better 



matching the spectral segments by shifting wavelength zero points in each order by < 0.05 
/im. For these pre-nominal Spitzer data, the SH wavelengths are more uncertain than relative 
intensity levels. We establish an absolute wavelength scale with the 17.03 /im 0-0 S(l) H 2 
line. We also improve the contrast of weak features by subtracting a continuum. 

We plot uncertainties in Figure [2] calculated from the root-mean-square difference be- 
tween the data and a seven-point moving boxcar average of the data. The statistical uncer- 
tainties, as measured from the difference between two individual integrations, are too small 
to be visible as error bars. 

Figure [2] clearly shows resolved emission features at 16.4, 17.4, and 17.8 /im. We detect 
weaker resolved emission features at 15.9 and 18.9 /im. The broadness of these emission 
features (0.1 - 0.6 /im, or 5 - 20 cm -1 ) implies a solid-state material, a molecular band, or a 
molecular vibrational mode broadened for the same reasons that the AEFs are broad. The 
only unresolved emission line in Figure [2] is the 0-0 S(l) H 2 line at 17.03 /im. 



3.2. Long-Slit Spectra 

We present our LL2 spectra (14 - 21 /im; R = 80 - 128) in NGC 7023 in Figures El 
HI El and [61 We plot uncertainties derived in the same way as for Figure [21 except that a 
five-point moving boxcar was used. 

Figure [3] illustrates our new results for slit A. These spectra demonstrate marked spectral 
variability from the PDR traced by the northwest filaments, dominated by the 16.4 /tm 
emission feature, the broad 17.0 /im feature, and the 0-0 S(l) H 2 line at 17.0 /im; to the 
regions between the PDR and the star, where the 16.4, 17.4, and broad 17.0 /im features are 
strongest; to the southeast, where the 18.9 /tm feature is most prominent and the 16.4 /im 
feature is extremely faint. 



Figure H] shows our results for slit B, from lPaper ll . which have been re-calibrated (~7% 
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change) and are illustrated here for comparison. The spectral features vary with spatial 
position in slit B much as they do in slit A. The 17.4 /im feature is most easily seen in 
spectra where both the 16.4 and 18.9 /im features are present. 

Figures [5] and [6] show our new data for slits C and D, respectively. Figures [5] and [6] reveal 
the 16.4 and broad 17.0 /im emission features and H2 emission in the northwest filaments; 
a strong 18.9 /im emission feature both northeast and southeast of HD 200775; and then a 
resurgence of the 16.4 /im feature and H 2 emission at the most southern positions. 

Many of the observed positions on the sky coincide to within ~1 pixel (~5") for slits 
C and D. We overplot the independently reduced spectra, at coincident spatial positions, 
on Figures [5] and [6j These independently reduced spectra at similar (within ±5") spatial 
positions are the best assessment we can provide of the quality of our data, because we are 
primarily limited by systematics rather than by signal-to- noise. 

We highlight some of our main results in Figure [71 This figure illustrates two pairs of 
spectra at similar projected distances (d — 24 - 25" or d = 49") from HD 200775, but in 
different directions (east or west). The continuum has been subtracted, by fitting a parabola 
to points well-separated from spectral emission features. 

Figure [7] demonstrates that the LL2 spectra east and west of the star have striking 
differences even at the same projected distance from the star, and at the same gas phase 
(neutral or molecular). The two western spectra are in the well-studied northwestern region 
of NGC 7023 containing narrow, bright filaments marking the PDR front there. The two 
eastern spectra are typical of the poorly-studied southeastern r egion of NGC 7023, where 



we recently discovered the 18.9 /im interstellar emission feature (IPaper ll ). 



The 16.4 /im emission feature has the highest intensity (in MJy sr 1 ) of any emission 
feature in the 15 - 20 /im spectra of NGC 7023 (Figs. El H U andEJ. The 16.4 /im feature 
emission, however, with its strong spatial correspondence with the northwest filaments in 
NGC 7023, has much lower intensity east and southeast of the star. Thus the 16.4 /im 
feature becomes weaker than the 18.9 /im feature in the two eastern spectra of Figure [71 

Figure [7] shows that the 18.9 /im emission feature is very prominent in the eastern 
spectra, especially when compared to its relative weakness in western spectra at similar d. 
This difference is observed at the southeast and northwest PDR fronts, as traced by H2 
emission (d = 49"), and within the neutral H I gas {d = 24 - 25"). 

Figure [7] shows that the 15.9 /im feature is easily seen at d — 49" NW, hinted at in the 
spectrum of d = 24" NW, but is not detected at d = 25" NE or d = 49" SE. We also conclude 
from Figure [7| that the 17.8 /im feature is observed only where the 16.4 /im emission feature 
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dominates the spectrum of NGC 7023: d = 24 - 49 " NW. The 17.4 /im feature is apparent in 
all spectra of Figured Its relative strength is difficult to estimate at this spectral resolution, 
due to strong blending with H 2 , the 17.8 /im feature, and probably at least one other feature 
(see below). 

At the spectral resolution of LL2 (R = 80 - 128), it is difficult to measure the relative 
strengths of individual emission features, particularly in the crowded 16 - 18 /im spectral 
complex. Still, there must be additional emission between the 16.4 /im emission feature 
and the 17.4 /im feature in order to explain the overall 16 - 18 /im emission complex. 



Smith et al.l (120061 ) fit this additional emission in LL2 spectra with a broad (1.1 /im FWHM) 



feature at 17.0 /im. We attempted a similar fit to this additional emission in NGC 7023, at 
the resolution of LL2 spectra, and were unable to find a unique fit. We examine this broad 
emission, underlying the 15.9, 16.4, 17.4, and 17.8 /im features, further in the next section 
using higher spectral resolution. 



3.3. Profile Fitting to Observed Spectra 



We model our SH spectrum, in order to qu antify the FWHM of each emission feature. 
We use a program called LINER (jPoggd 120041 ) to fit multiple blended Gaussians to the 
SH spectrum we ob serve. We illustrate o ur technique in Figure [SJ We also ran fits with 
Lorentzian profiles (IBoulanger et al.l 119981 1. which are similar to the Drude profiles used by 
Smith et al.l (120061 ); our conclusions do not depend on whether a Gaussian or Lorentian 



profile is used. 

First, we use a routine within LINER to interactively choose ranges of wavelengths for 
the continuum, as shown in Figure [8j LINER requires a minimum of two wavelength ranges, 
and a minimum of two wavelength points per wavelength range. The continuum for LINER 
is equal to a parabolic fit to the spectrum in these wavelength ranges. 

Second, we use another routine within LINER to make initial guesses at the height, 
width, and center of each Gaussian that is used to model each feature within the SH spec- 
trum. We do not fit profiles to features whose peak intensities are less t han 5% of th e 
continuum intensity, taking this as the limit for possible calibration artifacts (jUchidal 120051 ). 



Third, we use an interactive plotting tool within LINER to refine our guesses for the 
Gaussian parameters of each feature. This routine plots the individual Gaussians and the 
sum of the multi-Gaussian fit, resulting from our refined guesses, plotted over the data. 
This routine also plots the residuals between the data and the fit to highlight where the 
model needs adjustment. This routine allows us to modify our refined guesses, by adding 
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or subtracting features, or by changing a feature's central wavelength, peak intensity, or 
FWHM, in response to the quality of the model fit and its residuals. It is during this step 
that broader features, such as the 17.0 /im feature, are better defined. We iterate within 
this routine, modifying the refined guesses, until we have fit the observed spectrum by eye 
as well as we can. 

Fourth, we use a routine within LINER which uses standard numerical techniques to 
search for best fit parameters for each feature (modeled as a Gaussian) that minimizes the 
difference between the entire model fit and the observed spectrum, based on the refined 
guesses we provide. If the routine does not converge on a solution that is similar to our 
refined guesses, we repeat the process, modifying our initial and refined guesses, until we 
and the software concur on the fit. 

We show the resulting fit, using this technique, for our SH spectrum in Figure [SJ We 
illustrate the wavelength ranges for the continuum points, the derived continuum, the in- 
dividual Gaussians, the sum of the multi-Gaussian fit, and the residuals plotted together 
with the original spectrum. The residuals are another measure of the uncertainty in our SH 
spectrum. 

We measure the central wavelength and full width at half maximum (FWHM) of emis- 
sion features in our SH spectrum (R = 600). We list these derived properties in Table 1. All 
emission features are resolved except for H 2 . 

The wavelength calibration is uncertain for our SH spectrum, because it was acquired 
during pre-nominal operations. We therefore examine seven SH spectra <8" from Position 
B, from AOR key 3871232 (pid 28). These spectra, obtained from a spectral map with short 
integration times per position, are noisier than our SH spectrum and so are not reduced 
beyond the pipeline products or presented here. These spectra, however, are sufficient to 
measure the central wavelengths and FWHMs of all emission features except the 18.9 /j,m 
emission feature. For the 18.9 /im emission feature, we select the ten spectra with the highest 
signal-to-noise 18.9 /im emission feature, from the entire spectral map, to measure the central 
wavelength and FWHM of the 18.9 /tm emission feature. Nine of these spectra are 15 - 42" 
southwest of the star; one is 10" southeast of the star. 
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TABLE 1 

Observed Interstellar Emission Features at A/AA = 600 in NGC 7023 



K 


FWHM a 




(fim) 


(jim) 


References for Discovery 6 


15.9 ± 0.1 


0.3 - 0.6 


Sturm et al. 2000 


16.4 ± 0.1 


0.1 - 0.2 


Moutou et al. 2000; Van Kerckhoven et al. 2000; Sturm et 






al. 2000 


16.6 ± 0.1 


0.4 - 0.9 


This paper 


17.03 


0.021 


0-0 S(l) H 2 


17.2 ± 0.1 


0.2 - 0.8 


This paper 


17.4 ± 0.1 


0.1 - 0.2 


PaDer 1 


17.8 ± 0.1 


0.3 - 0.4 


Beintema et al. 1996 


18.9 ± 0.1 


0.2 - 0.3 


Paper 1 



(a) Central wavelength and FWHM for all features measured from SH spectra (AA = 0.021 
/am). The 0-0 S(l) H 2 line at 17.03 /im, used to fix the zero point of the wavelength 
calibration, is unresolved. 

(b) Reference for first detection of emission feature. 

The very broad 16 - 18 /im emission plateau underneath most of the features in Table 
1 is not well fit by a single Gaussian at ~17.0 /im, as found in the ISM previously at 



lower signa 
resolution (jSmith et al. 



-to-noise flBeintema et al. 



2004 



1996 



Smith et al 



Van Kerckhoven et al.l l2000l ) or lower spectral 
2006). Our best fit is with two broad features, 



one at 16.6 ±0.1 /im and another at at 17.2 ±0.1 /im. We find a best fit with two broad 
features both for Gaussian and Lorentzian fits, although the central wavelengths shift: 16.6 
and 17.2 /im for a Gaussian fit, or 16.8 and 17.1 /im for a Lorentzian fit. The 16.6 /im 
feature might be a long wavelength wing to the 16.4 /im emission feature rather than an 
independent emission feature, analogous to the long wavelength wings of the 6.2 or 11.3 /im 
AEFs. 



Discussion 



4.1. Laboratory comparison to NGC 7023 spectrum 



To prepare for the launch of ISO, iMoutou et al.l (119961 ) published the 14 - 40 /im labo- 
ratory spectra of neutral PAHs. They considered a number of PAH mixes, and successfully 
predicted the 16.4 /im feature in the ISM. They did not, however, predict the 18.9 /im emis- 
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sion feature. They also did not predict any of the other emission features observed in NGC 
7023 and elsewhere in the ISM, such as the 15.9, 17.4, or 17.8 /im features, or the broad 17.0 
/im feature which we fit best by two features at 16.6 and 17.2 /xm. 



Peeters et al.l (120041 ) construct the spectra of two laboratory mixes of PAHs, from the 
extensive PAH spectral database at NASA Ames Research Center. They compare these two 
mixes to ISO-SWS spectra of two different AEF sources, and find good fits. One source is 
dominated by the broad 17.0 /im feature, with no other strong emission features. The other 
source, CD —42 11721, is dominated by the 16.4 /im feature, as are many of our NGC 7023 
spectra. 



Peeters et al.l (120041 ) overplot our SH spectrum (15.1 - 19.5 /im; R = 600; see iPaper 1 
and Fig. |2J of NGC 7023 on the ISO-SWS spectrum of CD -42 11721, to show the spectral 
similarities. Both are dominated by the 16.4 /im feature. They then plot a spectrum of a 
PAH laboratory mix on top of their CD —42 11721 spectrum, to show the promising match 
in features. But they do not directly compare their spectrum of a PAH laboratory mix 
to our SH spectrum of NGC 7023. We show this comparison in Figure El with the PAH 
laboratory spectrum courtesy of Peeters and her collaborators. The 16.4 /im feature in the 
laboratory mix is an excellent match to the nebular 16.4 /im feature in NGC 7023. The 
weaker substructure in the laboratory mix bears some resemblance to the SH spectrum of 
NGC 7023. Figure [9] shows, however, significant mismatches when the laboratory and SH 
spectra are compared in detail. 



Peeters et al.l (120041 ) suggest that the 18.9 /im feature is due to ionized PAHs, based 
on fits to laboratory mixes of PAHs, but they do not specifically search for a PAH mix 
that would reproduce the spectrum of NGC 7023. Furthermore, they do not explore a 
mix including both neutral and ionized PAH species, as would be expected for NGC 7023. 
They plan a more detailed analysis of their laboratory data in the future. Their excellent 
match to the strong 16.4 /im emission feature, and their modest success in producing weaker 
emission features at 15 - 20 /mi, holds out promise that further laboratory work may provide 
identifications for some of the emission features at 15 - 20 /mi. 



4.2. Spatial variations in different emission features 

We observe marked changes in spectral features across NGC 7023 (Figs. El HI El El an< ^ 
[7]), which can be studied reliably using the long-slit capabilities of the IRS on Spitzer. In 
Figure [TUl we have divided the 16.4 and 18.9 /im feature intensity data into nebular regions 
east and west of HD 200775, and then plotted the spatial distribution of each. The statistical 
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uncertainties are less than the sizes of the points in Figure [10] the uncertainty in this plot 
comes from the intrinsic scatter between intensity values at different nebular positions with 
similar projected distances from the star. 

Figure [10] shows the 16.4 /im feature is strongly pe aked ~36" west of HP 200775. This 



is near the PDR front traced by H2-emitting filaments (ILemaire et al.lll996l ; lAn fc Sellgren 
20031 : IWitt et al.l[2006h northwest of the star. 



The 3.3, 6.2, and 11.3 /im AEF emiss i on is strongly peaked on t he northwest filaments 



( Cesarskv et al.lll996l : lAn k Selkren]l2003h . iMoutou et all fl2000h and lVan Kerckhoven et al. 



( I2000l ) identify the 16.4 /im feature with C-C-C vibrational modes of neutral PAHs, while 
the AEFs are widely attributed to PAHs. The spatial distribution of the 16.4 /im emission 
feature thus supports the current identification of it with PAHs. 

Figure [10] also shows that, by contrast, the 18.9 /im feature is symmetrically distributed 
east and west of HD 200775, with a smooth fall-off in intensity with d, the projected distance 
from HD 200775. We can fit a power-law to the distribution of the 18.9 /im feature intensity, 
I, with d. We find that I ~ d~ 15 . 

We find the same results for both the 16.4 and 18.9 /im features if the data are instead 
split into northern and southern sections: the 16.4 /im feature peaks in the north at the 
PDR front, and the 18.9 /im feature peaks symmetrically on the star. The feature intensities 
in Figure [JD] come from the area under a locally-defined linear baseline, rather than the 
Gaussian fit to the entire spectrum described in §3.31 but the results described in this section 
are unchanged if feature intensities from Gaussian fits to the entire spectra are analyzed 
instead. 

A basic difference between the 16.4 /tm feature and the 18.9 /im feature in NGC 7023 
is that the 16.4 /tm feature is bright in the northwestern PDR, where the column density is 
known to be high, while the 18.9 /im feature remains bright ~30" east and south of the star, 
where it appears the column density is low. 

The regions east and south of the star appear, at least in pr ojection, to lie inside the 
fossil bipolar outflow observed in the 24 /im Spitzer MIPS emission ( iPaper ll ). The northwest 
PDR, by contrast, appears to be part of the remnant molecular material surrounding the 
original bipolar outflow. It is possible that shock processing plays a role in the creation of 
the 18.9 /im feature carrier. 
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4.3. Identifying the 18.9 Micron Interstellar Emission Feature 

4-3.1. Minerals 



The strength of the 18.9 fim feature near the star, and its rapid falloff with d, leads 
us to consider whether the excitation of this feature could be equilibrium thermal emission 
from large grains, rather than non-equilibrium emission from tiny grains, PAHs, or other 
large molecules. We can (with difficulty) construct a model of equilibrium thermal emission 
that explains the 18.9 /j,m feature itself, but cannot construct a equilibrium thermal emission 
model for the underlying 18.9 fim continuum. In this case, the observed 18.9 fim continuum 
would be due to non-equilibrium emission, from tiny grains or large molecules. 

In the equilibrium thermal emission model for the 18.9 /xm feature, its carrier is large 
grains of unknown composition. We searched the literature for astrophysically important 
minerals, and also searched the Jena - St. Petersburg Database of O ptical Constants at 
|http : / / www, astro . uni- j ena.de /Laboratory / D at abase /j p do c/ index, html ( Henning et al.fl9~9~9 ) 
We can rule out th e following substance s as carriers of the 18.9 /iin feature: graphite 
(fortune &; Le3ll984h. a morpho us carbon (IColangeli et al.lll995f). hvdrogenate d amorphous 



carbon ( Ogmen fc Duley 1988 ). hydrocarbon nanoparticies ( Herlin et al. 1998h. amorphous 



silicates (both o 



Dorschner et al 



Koike et al 



ivines, Mg 2 , ; Fe 2 ,,. 



1995 



1993 



Koike et al. 



Jager et al. 



199, 



99c 



. 1999), 



1994 



S1O4, and pyroxenes, Mg^Fei-^SiOs; iJager et al.l 11994 



2000) , crystalline s i licate s (both olivines and pyroxenes 



'osch et al. 1999), corundum (a— A1 9 



Fabian et al 



Henning et al 



Koike et al.M2000L l2003h. qu artz (Si0 2 ; Koike & Shibai 



aaii. c ur unu u m \ a— aiqw-?; iPosch et alll999 ), spinel (MgAioW/i: i-ruscn et 
200 lh . FeO jHenning et all Il995l : If lenning fe Mutschke! Tl997l), MgO 



1995h. rutile fTiO ? :IPosch et ali ll999h. SiC JMutschke et al 



20011 ) . TiC (IHenning fc Mutschke 



FeCOs and CaFe(C0 3 ) 3 (IKemper et aL 



2001 



2002 ). crystalline or amorphous SISq (IBegemann et al 



,0 



Posch et al. 



19991 : iHenning fc Mutschke 
), car bonat es such as CaMgfCO.sW Ca CO.s, MgCO.s 



19971 ). or silicon nitride (Si 3 N 4 ; Idement et alll2005h . 



1996 



Kraus et al. 



19971 : 



2002h. FeS JHenning fc MutschkeHl997t). MgS JHonv et al 



Henning &: Mutschke 



4-3.2. Aromatic Molecules 

If the 18.9 /im feature is not due to emission from large grains in thermal equilibrium, 
then it, like the 16.4 /nn feature, is due to non- equilibrium emission from nanoparticies: tiny 
grains, PAHs, or other large molecules. Here we discuss one main class of large molecules: 
aromatics such as PAHs. 
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Rapacioli et al.l (120051 ) combine ISO-CAM circular variable filter (CVF) images of NGC 
7023 with laboratory and theoretical spectra of neutral and ionized PAHs, and analyze these 
with a decomposition technique, to construct what they call images of the ionized PAHs 
(PAH+) and neutral PAHs (PAH ) in NGC 7023. They measure the strength of PAH+ 
and PAH emission along a line from the star, through the northwest filaments, into the 
molecular cloud. Their results show that PAH + does not peak on the star, but rather peaks 
between the star and the northwest filaments, at d — ~24" from the star. They find that 
PAH peaks at d = ~42" from the star, in the northwest filaments. 



Flagey et al.l (120061 ) use ISO-CAM CVF images plus Spitzer IRAC images to derive 
improved images for PAH + and PAH in NGC 7 023, with a similar dec o mposition technique . 
Their images show the same results as those of iRapacioli et al.l (120051 ). iFlagey et al.l (120061 ) 
find that PAH + peak s midway between t he star and the northwest filaments, at d = ~22", 
not on the star itself. IFlagey et al.l (120061 ) observe PAH to peak on the northwest filaments, 



at d 



AO". 



We compare t he spatial distribution of our 18.9 /im emission feature intensity to the 
Flagey et al.l (120061 ) results in Figure [TT1 We calculate PA H + and PAH relative intensities, 
from images kindly provided to us by IFlagey et al.l (120061 ). at the same spatial positions at 
which we measure the 18.9 /im feature intensity. Figure [11] vividly illustrates how unique the 
18.9 /xm feature d istribution is c ompa red to the distribution of neutral and singly ionized 
PAHs inferred by IFlagey et al.l (120061 ) . Figure [TT] shows that the peak 18.9 /im feature 
emission is centered at d — 0, in contrast to the regions of peak PAH + (d = ~22" west) and 
PAH (d = ~40" west) emission. This is true even if one excludes the 18.9 /im data point 
at the position of HD 200775. 

Figure [TT1 clearly demonstrates that the distribution of intensity for the 18.9 /im feature 
we observe is spati ally quite distinct fro m the distributions for PAH + and PAH derived by 
Flagev et all (bood ). A similar plot of the IRapacioli et al.l (boosh PAH+ and PAH data yields 
the same result. PAH + and PAH are appealing carriers for the 18.9 /im feature, because 
of their likely ubiquity in space and the probability that laboratory data and theoretical 
modeling o f mixes of their spectra at 1 5 - 20 /im wil l be p ublished in the near future. If the 
analyses of IRapacioli et al.l (120051 ) and IFlagey et al.l (120061 ) are correct, however, concerning 
the location of PAH+ and PAH in NGC 7023, then PAH+ and PAH would both be unlikely 
carriers for the 18.9 /im feature. 



Witt et all (|2006|) observe the the Extended Red Emission (ERE) in NGC 7023 to be 



bright in the northwestern filaments. They show that the ERE is coincident in most cases 
with the H2 emission in the filaments. The 16.4 /im feature pea ks in these filaments (Fig. 
ITUl) . as do the AEFs (jCesarsky et al.lll996l ; lAn fc Sellgrenl 120031 ). This argues against the 
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ERE carrier and the 18.9 /im carrier being the same. IWitt et al.l (120061 ) propose that PAH 
dications (PAH ++ ) are the carriers of the ERE. This is a tentative argument against PAH ++ 
being the carriers for the 1 8.9 /im emission feature, but the PAH ++ identification model for 
the ERE (IWitt et al.l 120061 ) needs to be confirmed. 

If a mix of neutral and ionized PAHs does not provide a match to our NGC 7023 
spectra, then there are many other PAH variants to explore. A long list of references is 



given in lHalasinksi et al.l (120051 ). including (among others) PAHs with side groups, nitrogen- 



substituted PAHs, dehydrogenated PAHs and protonated PAHs. 



4.3.3. Fullerenes 



The 18.9 /im feature has a spatial distribution quite distinct from that of the 16.4 /im 
feature in NGC 7023. Figure [7| illustrates that spectra west of the star are dominated by 
the 16.4 /im feature (likely due to PAHs), while spectra east of the star are dominated by 
the 18.9 /im feature. Figure [TUl shows that the 18.9 /im feature peaks on the star, with an 
intensity proportional to d~ 1,5 ; the 16.4 /im feature shows no dependence on d. This spatial 
distinction between the 18.9 and 16.4 /im features suggests that two different species may 
be needed, rather than two ionization states of similar species. 



Jural (12004 ) has suggested that the 17.4 and 18.9 /im features in NGC 7023 are due t o 
neutral Cm- Cm w ould emit by non-equilibrium emission in NGC 7023 (IMoutou et al.lll999l ). 
Kwok et al.l (119991 ) detect narrow features at 17.85 and 18.90 /im in the C-rich protoplanetary 
nebula IRAS 07 135+1005. They tentatively detect a third narrow feature at 19.15 /im. 
Kwok et al.l (119991 ) no te that neutral Cm has four fundamental vibrational frequencies, at 7.1, 
8.6, 17.5 and 19.0 /im (IFrum et al.lll99ll ). They stop short of identifying the features observed 
in IRAS 07135+1005 as due to Ceo, however. The emission feature in IRAS 07 135+1005 at 
17.85 /im is not a good wavelength match for the Cgo 17.5 /im line. Furthermore. iKwok et al. 
(119991 ) do not detect the 8.6 /im C 60 feature in IRAS 07135+1005 (they are unable to search 
for the 7.1 /im line). 

We observe emission features at 17.4 and 18.9 /im, together, to the east and southeast 
of the star. These two features are a better wavelength match to the 17.5 and 19.0 /im lines 
of neutral C 60 . This suggests that these two features might be due to C 60 . The 17.4 /im 
feature, however, is also observed northwest of the star, where the 16.4 feature peaks and 
the 18.9 /im feature is weak. This casts doubt on attributing the 17.4 /im feature, at least, 
to C 60 . 



Two possibilities exist to explain why the 17.4 /im feature is observed in both east 
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and west of the star, while the 18.9 /mi feature is not. First, the 17.4 /im feature could 
be unrelated to the 18.9 /mi feature. In this case, neither feature is due to C 60 , and the 
17.4 /im feature has a spatial distribution which is distinct from both the 16.4 and the 
18.9 /im features. Second, the 17.4 /im feature could be a blend (at the lower resolution 
of our LL2 spectra) of two features. In this case, one contributor to the 17.4 /im intensity, 
perhaps Ceo, follows the 18.9 /tm feature, while the other contributor to the 17.4 /im intensity, 
probably PAHs, follows the 16.4 /im feature. These two possibilities might be distinguished 
by comparing higher resolution 15-20 /im spectra east and west of the star. 

Another way to test the possibility that the 17.4 and 18.9 /tm features are from Cm is 



to se ar ch for emission from its other two vibrational modes at 7.1 and 8.6 /tm (IFrum et al. 



199ll ). iMoutou et al.l (ll999l ) searched for the 7.1 and 8.6 /tm lines of C 60 at the nebular peak 
of NGC 7023, near the northwest PDR, and placed strong upper limits on C 60 . The 18.9 
/im feature is weak in this northwest region, however, so this does not rule out Ceo as an 
identification for the 18.9 /tm feature. By obtaining 7-9 fim spectra east and southeast 
of the star, where the 18.9 /tm feature is observed, one could search for the 7.1 and 8.6 /im 
lines of Ceo, to directly test the hypothesis that the 18.9 /im feature (and part of the 17.4 
/im feature intensity) is due to C 60 . 



Moutou et al.l (ll999l ) use the models of lBakes & Tielend (119941 ) to show th at fullerenes in 



NGC 7023 should be divided between neutral Ceo (35%) and ionized C^ (50%). IMoutou et al. 
( 119991 ) also put stringent upper limits on C^ at the NG C 7023 nebu l ar pe ak , from the ab 
sence of the 7.1 and 7.5 /im vibrational modes of Cg (IFulara et al.lll993l ). IMoutou et al. 



(119991 ). however, did not search for these emission lines east or southeast of the star. Spectra 
of NGC 7023 at 7 - 9 /im, at these positions where the 18.9 /im feature is bright, would 
simultaneously search for the 7.1 and 7.5 /im lines of Cg as well as the 7.1 and 8.6 /im lines 
of Ceo- No laboratory spectra at 15 - 20 /im exist for Cq Q , so it is not known where its longer 
wavelength lines might be. If Cgo is present in NGC 7023, however, then C^ should also 
be present at a similar abundance. Laboratory spectra at 10 - 40 /im of Cg are needed, to 
measure the wavelengths of all four infrared lines, and to measure their relative strengths. 



C70 has lines at 17.4, 18.1, and 19.0 /im ( INemes et al.lll994f). The lines at 18.1 and 
19.0 /im have equal intensity in gas-phase C70 emission ( INemes et al.l 1 19941 ). The absence of 
a detectable 18.1 /im feature in NGC 7023, particularly in the regions where the 18.9 /im 
feature is observed to be strong, makes C70 an unlikely identification for the 17.4 or 18.9 /im 
features. 
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4-3.4- Diamonds 



Many authors suggest that interstellar nanodiamond grains, similar to pre-solar nan- 
odiamonds extracted from meteorites, c an explain some infrared emission features seen in 



the ISM or circumstellar environments ( Koike et al.l Il995l ; iHill et al.l 119971 ; iHill et al.l 11998 



Guillois et al.lll999l ; IVan Kerckhoven et al.ll2002l ). Pre-solar meteoritic nanodiamonds have 
a me dian radius of 1.3 n m, and thus can radiate by non-equilibrium thermal emission in the 
ISM fcoike et al.lll995h . 



Laboratory spectra of pre-solar nanodiamonds have problems with surface contami- 
nation of the nanodiamonds, which adds strong features that can mask the weak fea- 
tures of diamond. Pure diamond itself has no infrared- active lines at 15 - 20 /an, 
but impurities and crystal defects break the symmetry and allow lines to appear in the 
diamond spectrum. Laboratory spectra of pre-solar nanodiamonds (from the Allende, 



(Lewis et al. 


19891: 


Co 


aneeli et al. 


1994 




Koike et al. 


19951: 


— , „ 
Mutschke et al. 


1995 


1997 




Andersen et al. 


1998: 


Braatz et al. 


2000l: 


Mutschke et al. 


2004 


) show no 



18.9 /im and thus rule out nanodiamonds as identifications for the 18.9 /im feature. 



Oomens et al.l (120061 ) argue that the identification of interstellar diamonds in meteorites 
make it worth considering smaller diamonds, such as diamondoids, as possible ISM com- 
ponents. Diamondoids, whose surfaces are fully hydrogenated, are the smallest diamond 
molecules. These diamond molecules can have a varying number of diamond cages, such as 
one (adamantane; Ci H 16 ), two (diamantane; Ci 4 H 2 o), or three (triamantane; Ci 8 H 2 4). The 
3-15 /im spectra o f diamondoids with o ne through six diamond cages have been measured 



i n the laboratory by lOomens et al.l (120061 ) . No laboratory data exist at 15-20 /im. iLu et al. 
(120051 ) predict vibrational frequencies at 15 - 20 /im for diamondoids with one through ten 
diamond cag es. They do not predict any features at 18.9 //m for the first five diamond cages. 
Piralil (120061 ) , who predicts vibrat ional frequencies at 15 - 20 /im for the same diamondoids 
measured by lOomens et al.l (120061 ) . predicts a 19.1 /im line for [121]tetramantane. This 19.1 
/im line is accompanied by a much stronger 16.3 /im line, however, which does not explain 
the strong 18.9 /im and faint 16.4 /im feature emission to the east and southeast of the star. 
Laboratory spectra of diamondoids at 15 - 20 /im are needed, but it looks unlikely that 
diamondoids can explain the 18.9 /im feature. 
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4.4. Further work 

Further observations, and further laboratory work, are needed to better identify the 
18.9 /im emission feature in NGC 7023, as well as other emission features at 15 - 20 /im. 

Further observations are needed where the 18.9 /im feature is bright in NGC 7023, east 
and southeast of the star, particularly high-resolution spectra at 7 - 9 /im and 15 - 20 /xm. 
Spectra should be obtained with enough spectral resolution to cleanly measure the strength 
of the 17.4 /im feature (LL2 spectra are too low resolution), and compared east and west 
of the star, to search for any blending of the 17.4 /im feature. This tests the hypothesis 
that the 18.9 /mi feature (and part of the 17.4 /im feature intensity) is due to Ceo- New 7 
- 9 /im spectra in regions where the 18.9 /im feature dominates the spectrum are needed to 
test the Cgo identification, by searching for the 7.1 and 8.6 /im lines of Cgo- This will also 
simultaneously search for the 7.1 and 7.5 /im lines of C^ . High spectral resolution (R ~ 
2000) is important at 7 - 9 /im because the C 60 features are superposed on the very strong 
7.7 /im AEF and are expected to be weak. 

The firm detection of t he 18.9 /im feature in only two out of a sample of 27 H II nuclei 



galaxies (ISmith et al.l 120061 ) demonstrates that this feature exists outside NGC 7023, but 
it is not common. Detailed observations of the phenomenology for all of the 15-20 /im 
features, but especially the 18.9 /im feature, are important to obtain in sources other than 
NGC 7023. 

Published 15 - 20 /im laboratory spectra of PAHs, fullerenes (especially Cg" ), and 
nanodiamonds/diamondoids are essential; these species, so far, are the leading candidates 
for the 18.9 /im feature. But the spatial dependence of the 18.9 /im feature intensity / in 
NGC 7023, with / ~ g?~ L5 (Fig. [10]), suggests that ionized and photochemically altered forms 
of these materials should be especially pursued. 

Alternatively, the 18.9 /im feature carrier could be a photodissociation or shock product, 
or a species only able to exist or survive near the star. Theory aimed at understanding the 
spatial distribution of the 18.9 /im feature intensity in NGC 7023 may be able to predict 
new ISM components for which 15 - 20 /im spectra laboratory spectra could be obtained, 
and compared with 15 - 20 /im spectra of NGC 7023. 

We are interested in identifying all the features at 15 - 20 /im, not merely the 18.9 
/im feature. In this context, laboratory 15 - 20 /im spectra should be obtained of PAH 
variants which have proven fruitful for studies of the AEFs, such as PAHs with side groups, 



nitro gen-substituted PAHs, protonated PAHs, or dehydrogenated PAHs (IHalasinksi et al. 



2005|). 
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A final possibility is that nanoparticles (radius ~1 nm) of minerals thought or known to 
be present in the ISM, such as silicates or SiC, might have spectra that differ from those of 
the solid-state mineral. Such nanoparticles would emit by non-equilibrium emission in NGC 
7023. Infrared gas-phase emission spectra of 1 - 2 nm radius nanoparticles, when technically 
feasible in the laboratory, would be an important area of exploration for ISM studies. 



5. Conclusions 



We observe the 15 - 20 jum spectra of the photodissociation region (PDR) and reflection 
nebula NGC 7023. We analyze four LL2 slit pointing positio ns in NG C 7023, including a 
re-reduction of the one LL2 slit pointing position presented in iPaper ll . We extract spectra 
at eight spatial positions alo ng each L L2 long-slit. We fit blends of Gaussian profiles to our 
observed SH spectrum from IPaper ll . to derive central wavelengths and FWHMs for each 
interstellar emission feature. 

We observe spectral emission features at 15.9, 16.4, 17.4, 17.8, and 18.9 /im. All of these 
features are resol ved in our SH spectrum . We observe the b road 17.0 Aim emission featur e 
first discovered by Beintema et al.l (1996), a nd later studied by I Van Kerckhoven et al.1 (120001 ). 



Smith et all (12004 ). and lSmith et al.l (120061 ) . We find, however, that it is better fit at higher 
spectral resolution by two features at 16.6 and 17.2 /im. We observe 0-0 S(l) H 2 emission, 
particularly at the filaments northwest of the star which mark the edge of the PDR there, 
and again at the PDR southwest of the star. 

The 16.4 /im feature is the brightest emission feature northwest of the star, and peaks 
in the northwest filaments. Southeast of the star, by contrast, the 18.9 /im feature is the 
brightest emission feature. The 18.9 /im feature peaks on the star, with its intensity falling 
off with projected distance d from the star as d~ 1,5 . The 17.4 £tm feature is present both 
northwest and southeast of the star. 



Moutou et al.l (120001 ) and IVan Kerckhoven et al.l (120001 ) identify the 16.4 /im feature 



with C-C— C vibrational transitions of neutral PAHs. We observe that the 16.4 /im emission 
feature peaks ~36" northwest of the star: this is in the PD R traced by H2-emitting filaments 



flLemaire et al.lll996l : lAn fc Selkrenll2003 



also peak in the northwest filaments (jCesarsky et al. 



Witt et al 



2006) . The 3.3, 6.2, and 11 .3 /im AEFs 



199ft lAn fc Sellgrenll2003h . The AEFs 



are widely attributed to PAHs, thus supporting the current laboratory identification of PAHs 
with the 16.4 /im feature. 



] 


3 eeters et al. 


(2004) 


(200, c 


')) and 


Flaeev et al. 



Peeters et al.l (120041) propo se the 18.9 /im feature is due to ionized PAHs. iRapacioli et al. 
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PAHs (PAH ) in NGC 7023. PAH+ peaks at d = 20 - 22", between t he star and the 
northwest filaments; PAH ° peaks at d — 40 - 42", on the northwest filaments (jRapacioli et al. 



20051 ; iFlagey et al.ll2006l ). We observe the 18.9 fxm feature to peak at d — 0", and to be 



symmetric ally distributed around the star with no peak at the northwest filaments. If the 
analyses of lRapacioli et al.1 (120051 ) and lFlagey et al.l (120061) are corr e ct, th en PAH + and PAH 
are unlikely as identifications for the 18.9 /im feature. IWitt et al.l (120061 ) propose that PAH 
dications (PAH ++ ) are the carriers of the Extended Red Emission (ERE) in NGC 7023, but 
they observe the ERE to be bright in the northwest filaments. If their identification of the 
ERE with PAH ++ is correct, then that would make PAH ++ unlikely as an identification for 
the 18.9 /im feature. Other PAH variants remain to be explored. 



Jural (120041 ) notes that neutral C 60 has features at 17.5 and 19.0 /im, while we observe 



emission features in NGC 7023 at 17.4 and 18.9 /im. Southeast of the star, we see these 
features together, with the 18.9 /im feature always stronger than the 17.4 /im feature. North- 
west of the star, however, we observe the 17.4 /im feature alone or much stronger than the 
18.9 /tm feature. Either the 17.4 and 18.9 /im features are not due to Ceo, or the 17.4 /im 
feature is due to Ceo plus a second 17.4 /tm carrier with a different spatial distribution. 



Moutou et al.l (Il999l ) suggest that ~50% of the C 60 in NGC 7023 should be singly ionized 



in NGC 7023. The wavelengths for the two longer wavelength infrared active C^" modes are 
not known; they are likely to fall near the wavelengths for Ceo (i.e. 17.5 and 19.0 /tm), but 
need to be measured in the laboratory. The 7-9 /tm lines of C 60 and C^ have been searched 
for in NGC 7023, but not yet southeast of the star where the 18.9 /tm feature is bright; these 
observations also need to be made. 

We search the literature for other identifications for the 18.9 /tm feature. We do not 
find a match among all other minerals that have been, to our knowledge, been proposed as 
components of the ISM and for which we were able to locate 15 - 20 /im laboratory spectra. In 
particular, we rule out amorphous silicates, crystalline silicates, amorphous carbon, graphite, 
SiC, corundum, and diamonds. 

We desire a better understanding of the 15.9, broad 17.0, 17.4, 17.8, and 18.9 /im 
emission features we observe in NGC 7023. This awaits new and better laboratory data 
on PAH variants, Qj~ , other fullerenes, nanodiamonds, diamondoids, and related species. 
Laboratory spectra, and ISM theory as a guide, are essential to sort through the possible 
identifications of these features. Detailed observations of these features in sources other 
than NGC 7023 will determine how widespread each feature is, and clarify whether the same 
phenomenology is observed elsewhere. 
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Fig. 1. — We illustrate how the IRS slit (10"5 x 168") at each pointing position (sli t A, 
B, C, or D) overlays different parts of NGC 7023 on a 2MASS Jskrutskie et all bonfih K s 
image. We extract LL2 spectra at eight spatial positions from each long-slit, with position 
1 furthest southeast and position 8 furthest northwest. Offsets from HD 200775 range from 
72" east 30" south (slit D, position 1) to 50" west 61" north (slit C, position 8). Slit B, 
position 3 is HD 200775. 
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Fig. 2. — Continuu m-subtracted SH spectrum (15.1 - 19.5 /im; R = 600) of Position B 
(0"W 29"N; IPaper lh in NGC 7023. We measure the intensity I v (MJy sr" 1 ) from the entire 
entrance slit (4'/7 x 11'.'3). We illustrate three overlapping SH orders, in contrasting colors. 
The wavelengths of individual orders are shifted by < 0.05 fim to match adjacent orders. 
The zero point for the wavelength scale is set by the 0-0 S(l) line at 17.03 /xm. We plot 
below the spectrum, as a measure of our uncertainty, the root-mean-square (rms) difference 
between the spectrum and a 7-point moving boxcar average of the spectrum. 
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Fig. 3.— The 15 - 21 /im LL2 spectra (R = 80 - 128) of slit A. We measure the intensity 
I v (MJy sr _1 ) with an extraction box of 10f!2 x 15'.'3. We label each spectrum with its offset 
(") from HD 200775. We mark the wavelength of the 0-0 S(l) H 2 line (dotted line). We 
plot ±l-rms uncertainties, from the difference between the spectrum and a 5-point moving 
boxcar average of the spectrum. 
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Fig. 4. — The 15 - 21 /zm LL2 spectra (R = 80 - 128) of slit B. We measure the inten sity 



I v (MJy sr" 1 ) with an extraction box of 1C/2 x 15'.'3. These spectra from iPaper ll 



arc 



re-calibrated (change of ~7%) for consistency and shown here for comparison. We label 
each spectrum with its offset (") from HD 200775. The spectrum of HD 200775 (offset 
0"east 0"north) is in units of F u (Jy). We mark the wavelength of the 0-0 S(l) H2 line 
(dotted line). We plot ±l-rms uncertainties, from the difference between the spectrum and 
a 5-point moving boxcar average of the spectrum. 
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Fig. 5. — The 15 - 21 fim LL2 spectra (R = 80 - 128) of slit C (magenta). We measure 
the intensity l v (MJy sr _1 ) with an extraction box of 10'.'2 x 15'.'3. Where there is spatial 
coincidence to within ~1 pixel (~5") between positions in slits C and D, we have overplotted 
the slit D spectrum (blue) on the slit C spectrum. We label each spectrum with its offset (") 
from HD 200775. We mark the wavelength of the 0-0 S(l) H2 line (dotted green line). We 
plot ±l-rms uncertainties, from the difference between the spectrum and a 5-point moving 
boxcar average of the spectrum. 
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Fig. 6.- The 15 - 21 //m LL2 spectra (R = 80 - 128) of slit D (blue), We measure 
the intensity l v (MJy sr _1 ) with an extraction box of 10'/2 x 15'.'3. Where there is spatial 
coincidence to within ~1 pixel (~5") between positions in slits C and D, we have overplotted 
the slit C spectrum (magenta) on the slit D spectrum. We label each spectrum with its offset 
(") from HD 200775. We mark the wavelength of the 0-0 S(l) H2 line (dotted green line). We 
plot ±l-rms uncertainties, from the difference between the spectrum and a 5-point moving 
boxcar average of the spectrum. 
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Wavelength (urn) 

Fig. 7. — Intensity I v (MJy sr -1 ), after continuum subtraction, vs. wavelength (/im), for 15 

- 20 fim LL2 spectra of NGC 7023. Note that the relative strengths of the 16.4, 17.4, and 
18.9 /im emission features (wavelengths marked by dotted lines) depend more on nebular 
position (east or west) than projected distance d from the star. North of the star (top two 
panels), at the nebular intensity peak, spectra show strong 16.4 /im emission. East of the 
star (bottom two panels), the 16.4 /im feature is much weaker, and the 18.9 jum feature 
appears strong by contrast. The 17.4 /im feature is bright both north and east of the star. 
The 0-0 S(l) H 2 line at 17.03 /im (wavelength marked by (dotted line) appears at d — 49". 
For spectra at d — 49", we plot ±l-rms uncertainties calculated from the difference between 
the spectrum and a 5-point moving boxcar average of the spectrum. For spectra at d = 24 

- 25", we plot ±l-er uncertainties from the average of two independent spectra obtained at 
the same position. 
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Fi g. 8.— Normalized SH spectrum (15.1 - 19.5 //m; R = 600) of Position B (0" W 29" 
N; iPaper ll ) in NGC 7023 (upper histogram) . We measure the intensity I u (MJy sr _1 ) from 



the entire entrance slit (4'/7 x ll'/3). We plot individual Gaussian profiles (solid curves) 
that have been fit to each feature (see text). We mark spectral regions used to define the 
continuum (vertical dotted lines), and illustrate the fitted parabolic continuum (solid line). 
We plot the final model spectrum (solid curve), which is a sum of the fitted continuum and 
the individual Gaussian profiles. We plot the difference between our SH spectrum and the 
model (lower histogram) as an estimate of our noise. 
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Fig . 9. — O ur continuum-subtracted SH spectrum of NGC 7023 (15.1 - 19.5 /zm; R = 600; 



sec 



Paper ll and Fig. [2] for details). We illustrate three overlapp ing SH orders, in con trasting 



colors. Overplotted (black curve) is a laboratory spectrum from lPeeters et al.l (12004 ). chosen 
to be a good fit to the ISO-SWS spectrum of CD -42 11721. NGC 7023 has an /SO-SWS 
spectrum similar to that of CD -42 11721 fjPeeters et alJbnO-d ). 
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Fig. 10. — The integrated feature intensities, /, for the 16.4 /im (top; filled squares) and 
18.9 /j,m (bottom; filled triangles) features, plotted vs. projected distance d (arcsec) from HD 
200775. We plot values to the east (positive d) and west (negative d) of the star separately. 
The 16.4 fxm feature intensity is strongly peaked in the west, at a value of d — ~36" west. 
The 18.9 jum feature intensity peaks on the star, at d — 0. We fit a power-law, / ~ <i~ L54 
(dotted curve), to the 18.9 /xm feature intensity. We measure the intensity (10~ 20 W m~ 2 
sr" 1 ) extracted in a 10"2 x 15"3 box from LL2 (R = 80 - 128) spectra. The angular 
resolution of each data point is 15'.'3. Statistical uncertainties are less than the point sizes; 
the uncertainties are due to scatter among different nebular locations at similar values of d. 
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Fig. 11. — Normalized intensity versus projected distance d (arcsec) f rom HP 200775, fo r 
the 18.9 /im emissio n feature (diamonds ; this paper), PAH + (plus signs; iFlagey et al.ll2006l ). 
and PAH (squares; IFlagey et al.l 120061 ) . We plot values to the east (positive d) and west 
(negative d) of the star separately. The 18.9 /xm emission feature peaks on HD 200775, at 
d = 0; PAH + peaks between the star and the northwest filaments, at d — ~22" west; PAH 
peaks on the northwest filaments, at d — ~40" west. The angular resolution of each data 
point is 15'.'3, 6", and 6" for 18.9 /im, PAH + , and PAH , respectively 



